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Abstract 



The conformal invariance of the low energy limit theory governing the electronic properties of 

graphene is explored. In particular, it is noted that the massless Dirac theory in point enjoys 

local Weyl symmetry, a very large symmetry. Exploiting this symmetry in the two spatial 

dimensions and in the associated three dimensional spacetime, we find the geometric constraints 

that correspond to specific shapes of the graphene sheet for which the electronic density of states 

is the same as that for planar graphene, provided the measurements are made in accordance to 

the inner reference frame of the electronic system. These results rely on the (surprising) general 

pH ' relativistic-like behavior of the graphene system arising from the combination of its well known 

special relativistic-like behavior with the less explored Weyl symmetry. Mathematical structures, 

^ such as the Virasoro algebra and the Liouville equation, naturally arise in this three-dimensional 

►^ , context and can be related to specific profiles of the graphene sheet. Speculations on possible 

^S| ' applications of three-dimensional gravity are also proposed. 

T-H ' 

o 

^ ; PACS No.: ll.30.-j, 04.62.-Fv, 72.80.Vp 

t^^ ' Keywords: Symmetry and conservation laws. Quantum fields in curved spacetime, Electronic trans- 

^^ ■ port in graphene 



1 Introduction 

In the low energy limit the electronic properties of graphene are well described by massless Dirac 
spinors ^ in two space and one time dimensions, hence an analog of a relativistic system but with 
characteristic velocity given by the Fermi velocity vp rather than the speed of light c [1] 

A = ihvF I d^x^j'^da^ , (1.1) 



where a = 1,2,3 is a (flat) spacetime index. This action is scale and conformally invariano 

The main messages of this paper are two, one general and one practical. The general suggestion 
is that intrinsic curvature of the two-dimensional graphene sheet within this special relativistic-like 
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^In the next Section we shall explain in detail the related concepts of scale, conformal and Weyl symmetries, both 
rigid and local, especially for the Dirac action. For the general case see [2]. 



behavior naturally leads to a general relativistic-like descriptiorO of ^I' regarded as a Dirac field in 
a curved three-dimensional spacetime back-ground 

A = ifivp I d^Xy/^^-f^Va^ , (1-2) 

where Va is the general covariant derivative, an important role being played by the local Weyl 
symmetry enjoyed by the electronic system. This suggestion, we shall see, comes very naturally 
when time is included in the analysis. The second, practical, main message of this paper is that the 
local Weyl symmetry of the graphene action can simplify the analysis of its electronic properties 
to the extent of allowing to find specific shapes of the graphene sheet for which the electronic local 
density of states (LDOS) is the same as that for a planar sheet, provided times and lengths are 
measured in accordance to the inner reference frame of the electronic system. 

We find that when the bending of the graphene sheet gives rise to a constant intrinsic two- 
dimensional curvature, the three-dimensional metric (the one that includes time) is truly curved 
also in the time direction, although, for the cases specifically considered here, only in a very special 
way (conformally fiat). Once such a general relativistic-like spacetime is accepted as physically 
meaningful, a special role is played by these conformally flat spacetimes 



i'J = V'(^)^i'J , (1-3) 



because in that case the curved space action (jl.2p is equal to the flat space action (jl.ip . provided 
the fields are related by a Weyl transformation 

A{^') = A{^) , (1.4) 

hence the Euler-Lagrange (Dirac) equations for ^' in the curved case are equal to those for ^ in the 
fiat case. As we shall see, this is so because of the special status of the massless Dirac action in any 
dimension. The three dimensions play a crucial role for other important matters, like for instance 
the particular form for the constraints for conformal fiatness (it is the Cotton tensor rather than the 
Weyl tensor that needs to vanish) but this particular instance, A^^') = A{^), is independent from 
the dimensions. This fact is known in field theory on curved spaces as an instance of conformal 
triviality [3J . We shall revisit that instance and shall show that it has the consequence of an invariant 
electronic LDOS. Our results do not seem to contradict existing studies of the nontrivial effects of 
curvature on the LDOS, such as, e.g, those reported in [4j, for two reasons: first, the ad hoc choice 
of the shape made there can well be outside the class of Weyl invariant choices (surely it is not 
within the class identified here); second, the measuring procedure adopted there seems to be a 
non-invariant one (we shall be more explicit later on this point). 

The paper is organized as follows. In Section 2 we shall recall the concepts of scale, conformal 
and Weyl symmetries in general and for the case of a massless Dirac field, and discuss the special 
case of conformal triviality and its impact on the Green functions. In Section 3 we shall relate 
the above to the case in point of the electronic properties of a graphene sheet and shall discuss 
the impact on the LDOS. In Section 4 we shall see which two-dimensional bending of the graphene 
sheet produces a case of conformal triviality in the three-dimensional spacetime and learn that these 
conditions are quite natural, i.e. constant intrinsic curvature of the two-dimensional sheet. This 
means three possibilities: i) the sheet is a plane; ii) the sheet is fiat but not planar (we find that this 
corresponds to two-dimensional conformal factors that are harmonic functions); iii) the bent sheet 



^The "like" sufSx here and in the following refers to the fact that this is not a fundamental theory but rather 
an effective one, one consequence being that the limiting speed is vf rather than c. Nonetheless, besides that, 
considerations are along the same lines as for the fundamental theories. 



is truly curved with a constant intrinsic curvature (we find that this corresponds to two-dimensional 
conformal factors that are solutions of the Liouville equation). There we notice that, although the 
system is truly a three-dimensional one, due to the nature of the metric used, two-dimensional 
conformal field theory structures, such as the Virasoro algebra or the Liouville equation naturally 
appear. Section 5 is devoted to collect our results and to discuss various possible directions for 
future investigations, in particular whether it makes sense for the physics of graphene to use the 
rich landscape of solutions of three-dimensional gravity, including black- holes. 

2 Scale, Conformal and Weyl Symmetries of the Dirac Action 

Scale, conformal, and Weyl symmetries are intimately related but different. For the general case of 
theories in any dimension and with arbitrary spin the issue was discussed at length in [2]. In this 
paper we want to focus on the special case of the massless Dirac field in two space and one time 
dimensions because it is relevant for graphene. To fully appreciate the peculiarities of this case, 
though, and for the sake of introducing notation and concepts let us start this Section by reviewing 
the results of [2j and by putting them in the right perspective for the applications we intend to 
pursue. 

Consider an action in n dimensions for the field <l>j, where i stands for a generic spin index, 
invariant under general coordinate transformations x^ — )• x^{x) (diffeomorphisms) and local Lorentz 
transformations Xa — )• (^a'^Xb 



Aig^,, <^i,Va<^i) = J dJ'x^Ci^,, Va^i) . (2.1) 



This action could be in flat space in curvilinear coordinates or in a truly curved geometry. We 
implicitly introduced the Vielbein e" and its inverse E^ 

riabelel = g^, e^E: = 5';, e^E^, = 5t , (2.2) 

the indices ^,1/ = 0, l,...,n — 1 that respond to diffeomorphisms (Einstein indices), while a,b = 
0, 1, ...,n — 1 respond to flat space transformations (Lorentz indices), rjab = diag(+l, —1, ...), y^ = 

det (7^;y| and the diffeomorphic covariant derivative that is 

V^^^ = ^^,<^^ + nf,/<l>j , (2.3) 

with Va = Ej^Vn, f]^/ = ^uj'^^{Jab)i , {Jab)i the appropriate form of the generators of the Lorentz 
transformationqj and 

u:,\ = el{5td, + Tl,)El, (2.4) 

is the spin connection obtained by requiring the full covariant derivative of the Vielbein to be zero 
(metricity condition) 

V^el = d^el - Tl.e'i + u:^ .e^ = , (2.5) 

where T^^ is the Christoffel connection. Torsion is defined as T^ny = T^r 1 = E\T"' where the 
expression in terms of Vielbein is T'^^^ = (5[^ ej,] + a;[^ Cj^])", see the Appendix. At this point it is 
not required that T^ = and to keep track of that might become important if topological defects 
described by a Cartan-Einstein model of gravity with torsion (dislocations) need to be included [5], 



^E.g., a $i — (j) is a, scalar field then Jab ~ 0, while if $i = 'i'a is a Dirac spinor, the case we shall be dealing with, 
then {Jab)i = |[7a,76]a . etc.. 



Due to diffeomorphism invariance of (|2.ip a rigid scaling of the coordinates can be transferred 
to the Vielbeiro (or, equivalently, to the metric). It is then clear then that when the flat space limit 
ofA{g) in ([21]), say 

A{r]ab, ^i, ^a^^) = J (i"x£($„ da<^i) , (2.7) 

is symmetric under rigid scaling 

x"" -^ e^x" and ^i -^ e'^^'^^i , (2.8) 

where d^ is the scale dimension of the field dictated by the kinetic term, then A{g) is symmetric 
under 

< -^ (i^e-l and $i -^ e'^'^'^^i . (2.9) 

The transformations ()2.9p are called rigid (being a constant) Weyl transformations. They are called 
local Weyl transformations if (t{x). 

Let us now recall that often Poincare and scale invariant actions are also invariant under the 
full conformal group of transformations, namely, on top of the Poincare and scale transformations, 
they are invariant under those transformations for which the scaling parameter iqj 

(^^c[x) = - ln(l + 2cx + c^x^) , (2.10) 

where Ca is a constant vector parameter, and SC stands for special conformal transformations. This 
is the case of the massless Dirac action of our concern here but it is not always true. 

The exact criterion to establish when it is true that scale invariance of the flat space action 
implies full conformal invariance of the latter invokes the local Weyl invariance of the curvilinear 
action, as one might suspect by noticing that ctsc is coordinate dependent. The general procedure 
is presented in detail in [2] and we shall not do it again here but it is crucial for us to understand 
the special status of the Dirac action in any dimension with respect to fields of integer spin. 

Thus, let us consider the case of a general x-dependent scaling parameter (t(x), i.e. not con- 
strained to be o"sc- For cr{x) the Weyl transformations (j2.9p are not symmetries of A{g) in ()2.ip . 
just like the scale transformations ()2.8p when a{x), in general, are no longer symmetries of A{r]) in 
(|2.7p . To make A{g) invariant we need to gauge it d la Weyl (these considerations now presume 
torsionless connections) 

A{g^,y,^iy a^i) ^ Aw{g^,v.^i,'Da^i) , (2.11) 



where 

and 

with 



V^<^i = V^^i + (Ar)/Ty,$,- (2.12) 



(A^,)/ = d^g^,5i + (J^,)/ , (2.13) 

W^^W^- G^ and V^^, ^ e'^^'P^^i , (2.14) 



Under diffeomorphisms 



dx 



\ 



~^l^i) = 9^6^(2;) , (2.6) 

which, requiring i^ = e~'^a;^, gives ej^(i) = e" e'^^^(x), that is the same result obtained with i^ = a;^ and requiring 
e",(a;) = e''eji(a;). The second choice, though, hides an important difTerence, namely that x in ej^ is indeed just x, 
hence the transformation now does not involve the coordinates but only the fields ej^. 



Other choices often used are ctsc = —\ ln(H-2ca:: + c^2;^) and ctsc = — ln(l + 2ca; + c^a::^) (and sometimes Ca 



-Ca 



as in [5]). They are all equivalent to the application to Xa of an inversion followed by a translation by Ca, followed by 
an inversion. 



under (|2.9p with a{x). Here u^j = d^a and J^,y = e" JabC^. As for usual gauge theories, the action 
is modified to a new one that is invariant: Awio) — ^ •^wid) under (j2.9p with cr(x). If and only if 
this produces an Awid) where W^ only appears in a combination, say it r2^^[WA]! such that in the 
flat space limit r2tt;,[(Tc] = has o"sc in (|2.10p as the only solution, then the flat space action (j2.7p is 
conformal invariant. In general this only happens for certain combinations of the spin of the fields 
and of the dimension of the spacetime. 

What just said needs only a little adjustment in relation to the two-dimensional case, where the 
o"scs are many more than just (I2.10p as all harmonic a{x)s (i.e. the infinite solutions of dad°'a = 0) 
correspond to conformal transformations. The procedure in this case gives as necessary and sufficient 
condition for a scale symmetric A{rj) to be fully conformal invariant that W^ enters Aw{g) as ^^. 
Special attention is needed for the scalar fields, as pointed out in [7J. 

In this respect, the status of the massless Dirac actioroin any dimension is very special becausqll 

Awig^.u:'^ ,Va^) = i j d^x^ ^^^EliiV ^ + K;:W,)^ = A{g^,,^ y a"^) , (2.15) 

due to 

7^A^^ = d^7^ + ^^'J^^ = , (2.16) 

where we used d^ = (1 — n)/2 and the definition of the Lorentz generators J^y = \['^^,'~iv\- This 
means that not only the Dirac action in flat space 

A{riab.^,da^)=i j d^X^^da^ , (2.17) 

is conformal invariant but that the curvilinear or truly curved space action 

A{g^,,, ^, V,^) = i J d^x^ ^rEJiid^ + ^oj^'Jbc)^ , (2.18) 

is local Weyl invariant, as can be also seen directly by transforming (j2.18p under (j2.9p with a{x) 
obtaining 

A{gf,„ *, V„^) ^ A{g^,, ^, V,^) + i J d^'x^ %'^A;a,^ = Aig^,, ^, V^^) . (2.19) 

Here it was used u^^ — )• w^"^ + (e^e[^ — e|^e"). This and other relevant formulae are in the Appendix. 

In any dimension this is a very large symmetry that A{g^u, ^, V^^) enjoys because there are no 
restrictions to cr{x) of any sort. This happens for A{g) as it stands which means that the theory is, 
we may say, "intrinsically gauged": Aw{g) = A{g). We shall exploit this symmetry in what follows 
for the graphene action where a crucial role is played by the dimensionality of the problem, which 
is two space and one time dimensions. Let us see here what this symmetry means in physical terms 
in any dimensions, including the case of graphene. 

When two metrics are related as g' = e^"^^' g^y the classical physics for the field ^' = e *'^*^^^^ 
in g' is precisely the same as the classical physics for the field ^ in g^y because A{gij,u, ^, Va^) = 
A{g' ^y,"^' jVa"^')- Furthermore, we can choose 

g'nu = V^.u , (2.20) 



''Although for shghtly different reasons, what foUows is also true for scale invariant actions of fields of any half- 
integer spin (Rarita-Schwinger). 

^For the moment, we shall take h = 1 = v with v a parameter with the dimensions of velocity, for graphene it is 
V = vf, the Fermi velocity, for truly relativistic Dirac it is, of course, v = c. 



and obtain 



I rx^^^'^EiiV^^ = ijdrx e-("-i)- '^7«(5a - ^^a)* = i j (f^x^'^^'da^' , (2.21) 



when 

5^. = e-2-(-)r?^. and ^ = e^-(-)^' , (2.22) 

where we used the simple algebra illustrated in the Appendix to write the second form of the action. 

We are then dealing with a conformally invariant field in a conformally flat spacetime, a case 
sometimes referred to as conformal trivialitv^ [3j. This means that the effects of curvature are null 
on the classical physics of a massless Dirac field ^ as in (I2.22P if the spacetime is only curved in 
a conformally flat fashion: if the metrics can only be the conformally flat we know from here that 
the classical physics is once and for all governed by the flat space action. Of course this is not 
always true as there is no reason in general to consider only conformally flat metrics, but when this 
happens there is the "magic" consequence that although the spacetime is curved this has no effects 
whatsoever on the classical physics of the system, provided the field is properly redefined. 

What for the quantum case? If we define the flat propagator as usual 

iS' {xi,X2) = '(0|vl/'(xi)^'(x2)|0)' , (2.23) 

where 

i,^,,5'(xi,X2) = 5"(xi,X2), (2.24) 

we can see the effect of the Weyl transformations by considering the unitary operator that imple- 
ments the transformations quantum mechanically, i.e. 

^(x) = U^'{x)U~^ = e^''(^)^'(x) . (2.25) 

By writing U = e and by using [^'(x), ^ (y)]+ = 6^{x — y), and from 

U^'ix)U~^ = ^'{x) + [B,¥{x)] + ^[B,[B,^'{x)]] + --- 

77—1 1/77 — 1 \^ 

= *'(x) + ^- aix)^'ix) + - [-^ a(x) j ^'{x) + • • • , (2.26) 

we easily obtain 

u = ^^v\^-^ j d^y<y)^'\yW{y)] , (2.27) 

i.e., 5t = -B and U{-a) = U~^{a). Notice also that [^,7"] = 0, hence U'j'^U-'^ = -f" and 
^(x) = U^'{x)U^^ = e^~'^'^-'^'(x). The vacuum associated to the Weyl transformed fields is of 
course 

|0) = C/|0)' = e^^il^ j d-ya{y)^'\yW{y)^ |0)' , (2.28) 

that shows the typical condensate structure of a quantum field in a curved (classical) background 
[8]. Now, as discussed at length, e.g., in [8], we can do two things: either we consider 

i5(xi,X2) = (0|^(xi)f'(x2)|0) (2.29) 

= ' {{)\U'^U-^' {xi)U-^UW {X2)U-^U\(S)' (2.30) 

= '(0|'I/'(xi)^'(x2)|0)' = i5'(xi,X2), (2.31) 



Notice, once more, that for fields of half-integer spin conformal triviality needs not a coupling of the fields to 
curvature tensors to have local Weyl invariance, while this is necessary for integer spin. 



or we consider 

i5"(xi,X2)= '(0|*(xi)^(x2)|0)' = e'^("(^i)+"(^2))/^o|^'(xi)^'(x2)|0)' (2.32) 

iS-%xi,X2)^ (0|'f'(xi)^'(x2)|0) = e^("(^'i)+"(^2)) '^o|^'(xi)^'(x2)|0)' , (2.33) 

notice that (0|^'(xi)^'(x2)|0) -> '(0|^(xi)^(x2)|0)' when a -^ -a, as it should be from U{-a) = 
U^^{a). Summarizing: 

S{xi,X2) = S'ixi,X2) (2.34) 

5'^(xi,X2) = e'^('^(^i)+'^(^2))^/(^^^^2) , (2.35) 

with S"^(xi,X2) the Green function for the curved Dirac operator 

i Wx^S^x^,X2) = ^r(xi,X2) , (2.36) 

as proved in the Appendix. 

The physical meaning of the two procedures is clear. In the first case, Eqs. (|2.29|) - (|2.3ip . the 
measurements are made within the frame of reference that is co-moving with the particles described 
by the field ^, hence the effects of curvature are completely removed also at the quantum level. In 
the second case, Eqs. ()2.32p and ()2.33p . the measurements are made by an observer that sees the 
(quantum) effects of curvature in the form of a condensate in the vacuum, Eq. (I2.28p . 

We faced this sort of situation in [8j where we dealt with vacuum expectations of the kind 
(0(e)|Ar|0(e)) (or viceversa (0|iV(e)|0)), with U{e)NU-^{e) = N{e) and U{e)\0) = |0(e)), where 
1 0(e)) is the vacuum for a scalar field, e.g., in a Schwarzschild spacetime, with e a function related 
to the acceleration a or to the surface gravity k (near the horizon the space-time is always Rindler, 
see, e.g., [9]). The number operator A^ is that of the same field in an inertial frame, i.e. a frame that 
is not freely falling into the black-hole with acceleration a. In other words, N represents a particles' 
"counter" placed asymptotically far-away from the horizon in a region that can be approximated 
as flat. It is this counter that sees the condensate structure of |0(e)) 

(0(e)|iV|0(6)) = sinh2 e{p) = ^^^ , (2.37) 

where p = (ri, k) is the four momentum in the curved frame, /? = 1/T is the inverse temperature, 
the sinh^ e{p) descends from the algebraic structure of U{e) and the relation of the latter to the 
Bose distribution comes from the minimization of the free energy [8]. If, on the other hand, we let 
our particle counter fall into the black hole, this will not see any horizon, nor there will be a way for 
it to detect the acceleration, as the principle of equivalence dictates. Mathematically this translates 
to 

(0(6)|A^(e)|0(6))=0, (2.38) 

that is precisely what we get in the absence of any gravitational/acceleration effect, (0|A^|0) = 0. 

Returning to our case, if we think of ^^ as our "counter", we can apply the same logic just 
illustrated, i.e. S{xi,X2) = S'{xi,X2) means that we see no effects of curvature, while S'^{xi,X2) 
contains the information on the quantum vacuum condensate. Two instances are important here. 
First, the physical picture can only become fully clear when the particular (conformally flat) metric 
9ij.u = e~^'^?7^;y has been specified. When we can do that we can give meaning to lengths and times 
measured according to the line element ds'^ = e~'^'^{dt^ — dx'^), hence we can appreciate the meaning 
of the invariance S{xi,X2) = S'{xi,X2)- What this means practically (i.e. what needs to be done 

7 



to the measuring apparatus) depends upon a. Second, in [^ we have been deahng with a truly field 
theoretical setting, i.e. with an infinite number of degrees of freedom and the related infinite volume 
limit. In that case the Bogoliubov operators C/(e) were not well defined in the limit and unitarity 
was lost (on this see, e.g., [10]). For the applications to graphene we can safely take the view that 
unitarity of U{a^ is never lost because the size of the sample is always finite and the number of 
degrees of freedom cannot be infinite because to excite more and more degrees of freedom we would 
need higher and higher energy than that necessary for the "Dirac-like" approximation to hold (see 
next Section). 

A last remark on the quantum regime is related to the fact that what just said seems to apply 
to any [/, i.e. not just to the C/((t)s in ()2.27p that generate the Weyl transformations hence are 
symmetries of the classical theory. To select which C/s are to be used one can look at the path 
integral version of the Green functions (see, e.g., [11]) and use the classical Weyl-symmetry of the 
action A{^,a) = A{^') 

PI . _ /PM/P'l'^(xi)^(a:2) exp{a(M/,a)} 

^ (^^'"^) = — jv^vw^M^M^:^)) — ('•'') 

^ zi^(.(.,)+.fe)) /P^^P'l'^ ^\x^)^\x^) exp{a(^I/0} 

/P^'P^'exp{a(^')} ' 

= e'^("(^i)+"(^2))S'(xi,X2). (2.41) 

Here, since (j{x) is an external field (see later), we used P^ = e^~'^^^-'P\E'', similarly for 2?^, hence, 
in the ratio the overall factor e~2~('^(^)+'^(s^)) cancels. From the above we see that the Green function 
computed via standard path integral methods, S*^^, is indeed S"^ in p. 320 . This is so because we 
can as well read (|2.40p as 

\V^'V^>' ^{x{)^{x2) exp{a(^0} _ 

7^^7^f7^^^^{-4(^7)} = (0|*(xi)^(x2)|0) . (2.42) 

Nonetheless, after having exploited the classical symmetries in the path integral Green functions we 
can then use the operator approach to find the C/s corresponding to these symmetries, for us here 
the U{a)s in (12.270 . and as well construct the invariant Green function as in (I2.30p . 

Let us now compare Weyl gauge symmetry to standard gauge symmetry. To this end consider 
the Dirac action 

A{T]ab,^,Aa)=iJrx^^^{da + Aa)^ (2.43) 

invariant under 

^ ^ § = e*"^^)'^ , ^ ^ ^ = e-^"(^)^ and Aa ^ Aa = Aa - iua , (2.44) 

where Aa is the field that needs to be introduced to adsorb the derivative of the gauge parameter 
daO = aa- The invariance means that A{r]ab, ^j ^a) = A{r]ab, ^, Aa), i.e. we can use the fields ^ and 
Aa or the gauge transformed ones ^ and Aa and we would not be able to see any difference because 
the physical effects of a nonzero a are null for this form of the action. For this part the analogy with 
the Weyl symmetry we just discussed is total: we can use the fields ^ and a or the Weyl transformed 
ones ^' and a' and we would not be able to see any difference. Furthermore, for the standard gauge 
theory (j2.43p there is no physical meaning we can ascribe to a: it is just the parameter of a kind of 
rotation of the field point by point. This is the reason why, besides topological non-triviality, the 
gauge field itself is not observable. Here the analogy with the Weyl case breaks down twice, once 
because a is directly related to the metric hence indeed it has a physical meaning on its own, and 



once because, for the case in point of an intrinsically Weyl-gauged action, the Weyl-gauge field is 
absent hence it is not observable because it is not there in the first place. Thus we can indeed take 
the view that the non-transformed field ^ is performing a curved space physics in a metric e~'^'^r], 
but identical physical results are obtained by considering the Weyl-transformed field ^' performing 
a flat space physics. 

Weyl symmetry is an internal symmetry hence, as such, the associated Noether current is of the 
form Ja = Hl^S^^i, where, as usual, IT* = 6A/Sda^i- Now we have to decide how to treat a in 
(I2.2ip . namely if we need to include it in the current as one of the fields <l>j or not. To properly treat 
o" as a dynamical field we would need to include a kinetic term for it, just like we do for standard 
gauge theories where we add to the gauged action a further term, usually of the form FabF , where 
Fab is the field strength. But the Dirac action in point needs not be Weyl-gauged, that is why 
there is no VF^ field, hence the associated kinetic term, if any, can only be relative to the geometric 
quantity a, hence to the metric. 

Thus, if no dynamical meaning can be ascribed to cr{x), i.e. cr(x) is an external field, the Weyl 
current is 

Ja = ie(i-")"'l'7a^ = i^'-fa^' . (2.45) 

The action (j2.2ip is also invariant under standard rigid gauge transformations, i.e. (|2.44p with 
Ua = 0, whose associated Noether current is 

ja = e^i-^^'^^Ta^ = ^'7a^' = -iJa • (2.46) 

The way we interpret this result is that we have to choose which current to retain as physicao, and 
we may as well choose ja- With this choice no further conserved current is introduced when rigid 
Weyl symmetry is present. The action is not invariant under local gauge transformation, a{x) but 
it is under local Weyl and the physical meaning of this symmetry we have explained earlier. 

It might look surprising that the Weyl current Ja has nothing to do with the energy-momentum 
tensor but this is as it should be. The Weyl transformations are abelian internal transformations 
hence no T^;^ (a current associated to the non-abelian spatiotemporal group SO{n,2)) could ever 
appear this way. The way to obtain the latter is through 

9 /\ A 

T'"' = — = i^^^V^)^ - i5^^^7^VA'f , (2.47) 

^/gogijiv 

whose trace is zero on-shell {JJ'^ = 0) with no needs to add improvement terms. The absence of 
such improvement terms is the effect on T^^^ we have to expect from local Weyl symmetry. 

3 Weyl-Gauge Symmetry of Graphene 

Graphene is a two-dimensional honeycomb lattice of carbon atoms arranged in two triangular sub- 
lattices, say La and Lb, whose electrons in the vr-bonds belonging to one sublattice can hop to the 
nearest neighbor sites of the other sublattice. The electronic properties of graphene are ascribed to 
these electrons. The elastic properties, instead, are ascribed to the a-bonds and involve an energy 
orders of magnitude stronger than that relative to the vr-bonds. Let us focus on the electronic 
properties which, in the low-energy approximation, are described by the Hamiltonian (we set h = 1) 

3 
H = -tY,Y. {a){r)h{f+ Si) + h\r + Si)a{r)) , (3.1) 



®One could also take a less conservative view and consider complex currents J'a — \J'a\e and ja ~ \ja\e^'^, with 
\J^a\~ \ja\ and X — ''r/2 = (j), but we shall not do it here. 




• = sublattice La o = sublattice Lb 

Figure 1: The honeycomb graphene lattice. 

where t ~ 2.7 eV is the nearest neighbors hopping parameter [12] (the next-to-nearest neighbors 
hopping parameter t' is taken to be zero, and that is the low-energy approximation), a (a^) and b 
(6' ) are anti-commuting annihilation (creation) operators for an electron in the La sub-lattice and 
in the Lb sub-lattice, respectively, and all vectors are two-dimensional r = {x,y) and Si = (sf,sf), 
as described in Fig. 1. The relevant anticommutation relations are 

[a{f),a^r')]+ = 5\f- P) , [b{f+ s,),b\P + Si)]+ = 5\f- P) , [a,b]+ = [a,b^]+ = (3.2) 



First one notices that due to the geometry of the hexagonal lattice the basis vectors of the 

sublattice La are 

P - P - 

(3.3) 



ai = -(\/3, 3) and 02 = -(\/3, -3) , 
and the vectors moving from a site of La to the three nearest neighbors Lb sites are 



si = £(0, -1), S2 = -(\/3, 1) and S3 = ^(-Vs, 1) . 



(3.4) 



In the following we shall set the lattice spacing i to 1. Taking the Fourier transform, a{f) 
'}2,ra{k)e'^^'^ , etc, the Hamiltonian (j3.ip becomes 



H = Y,{f{k)a\k)b{k)+\..c.), 

k 



where 



f{k) 






,^/3, 



-te-'^'y I 1 + 2e*2''" cos(^^A:^) 



The 1-particle spectrum being given by 



E{k) = ±\f{k)\ 



(3.5) 



(3.6) 



(3.7) 



the modes with zero 1-particle energy are easily found as solutions of f{k) = 0, i.e. k^ = (±—^,0), 
where the superscript D stands for Dirac (points). If we linearize around k^, k± ~ k^ + p, and 
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f±{p} = f{k±) = ±y{Px i ipy), a±{p) = a{k±), b±{p) = b{k±), then H in (j3.ip can be written as 



t;p5:[(at+,6l,) 



Px+ ipy 



P 

t J 



{bl,a 



Px - Wy 
Px+ iPy 







64 



Pa' - iPy 







m 



where vp = ^ is the Fermi velocity (at i = 1) that we shah now set to 1. If we define 



(3.8) 



^+^ 



and V- 



6_ 
a_ 



(3.9) 



and a = (ui, C72), with 



0-1 



1 

1 



0"2 



-i 
i 



(3.10) 



the usual Pauli matrices, and we Fourier-transform back to configuration space, then the Hamilto- 
nian (|3.8p becomes 

H = J2{'^+^-P'^+-'^-^-P'^-) =-i [ d'^x(ijl_a-d'il;+-i;la-dij-^ . (3.11) 

p 

The spin 1/2 structure here is entirely due to the lattice, that is why it is often called pseudospin, 
but we shall treat it as aproper spin, as suggested, e.g., in [13]. With the introduction of the four 
component Dirac spinoic^ 



^ 






of the following form for the Dirac matrices 



a 



a' 
-a' 



, /3 



6_ 


V "- 


-J 





h 




z = l,2. 



(3.12) 



(3.13) 



which lead to the usual definition of 7-matrices 

satisfying 

[7^7']+ = 2r?'^^/4 



(3.14) 



(3.15) 



with a, 6 = 0, 1, 2, the "fi? relativistic" Dirac Hamiltonian in (13. lip can be further compacted to 

H = f (fx(-i^j-d^) = f (fxn , (3.16) 



We prefer to work with four components Dirac spinors because we consider the graphene system as living in a 2 
+ 1 dimensional spacetime and Weyl spinors, although natural in the two (spatial) dimensions we have considered 
till now, do not exist in three (spatiotemporal) dimensions. 
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where * = ^t^o 

Thus, starting from the instance that the energy dispersion relations (|3.7p become hnear at 
the Dirac points, taking into account only the low energy contributions, including the effect of the 
honeycomb structure of Fig. 1 and considering situations that do not include defects or impurities, 
the effective description of the electronic properties of graphene is given by a special relativistic-like 
Hamiltonian from which through customary Legendre transformation, L = W^ "^ — Ti, a Lagrangian 
density can be derived. What is obtained is precisely the scale- and conformal-symmetric massless 
Dirac action in flat space we discussed in the previous Section with n = 3 

A = i f (fx^-i^da-^ . (3.17) 



Now, let us suppose that the graphene sheet is hent. By this we mean that the two-dimensional 
sheet of graphene is not planar, hence curvilinear coordinates are better suited. A typical situation 
we have in mind is, of course, that of the ripples [Hj, [15], [IB] (see also the review |17j), that are 
found in experiments and that are the subject of much theoretical and experimental investigations 
on their origins and their impact on the electronic properties of graphene, but other geometries are 
also within the reach of the following analysis. The graphene action to consider is then 



A = ijd\^ ^l^E^id^ + ^Lo^'Jbc)^ , (3.18) 



where g^u is a three-dimensional metric that bears the information about the bending in the two 
spatial dimensions (we shall be more precise about this in the next Section). Thus the graphene 
action, in the limits recalled above, is a particular instance of the "intrinsically Weyl-gauged" action 
for massless Dirac fields we have discussed before. 

Every bending of the two-dimensional sheet produces a conformally flat two-dimensional metric 
gap, Oi,l3 = 1,2 for the simple reason that every two-dimensional metric is conformally flat. It 
is matter of finding the coordinates system where this is evident, but for a well known general 
argument this is always possible. These conformally flat g'a/js describe graphene sheets of which: 
(i) one is planar, hence intrinsically (and extrinsically) flat; (ii) some are non-planar but are in 
fact intrinsically flat; (iii) some have constant intrinsic curvature; (iv) some have non-constant 
intrinsic curvature. Due to the local Weyl invariance in any dimension that we discussed in the 
previous Section, if the dimensionality of the problem was indeed two we would have immediately 
the result that in each one of the cases (i)-(iv) the electronic properties of graphene are described 
by a conformally trivial system and we could take advantage from this every time. 

On the other hand, the correct dimensionality is three because there is time. In the next Section 
we shall see that the presence of time, the special relativistic-like form of the action and the two- 
dimensional conformal flatness of the graphene sheet naturally lead to a general relativistic-like 
metric, namely a 2+1 dimensional conformally flat metric, including the truly curved. We resort 
to such a discovery while asking the question: is there a particular shape of the two-dimensional 
sheet that gives raise to a conformally trivial Dirac system also in three dimensions, i.e. including 
time? We postpone the answer to the next Section as we want to come back now to the effect of 
local Weyl invariance on the density of states of graphene. 

For graphene, the quantum version of the local Weyl invariance discussed in general in the 
previous Section can have striking consequences. First we notice that the electronic LDOS can be 
written in terms of the two-point function [18] (see also [19] and [1]) 

p{E, x) = -ImTi {g{E, f , f)7°) , (3.19) 
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where G stands for a generic Green function (S, S' , S'^ or S~'^ of the previous Section) and to obtain 
its energy dependence it is customary to move to the Hamiltonian form of the Dirac propagator, 
formahy: G{E) = 1/{H — E). Then we need to apply the general discussion to the case of graphene, 
by first asking whether it makes sense to have a truly 2+1 dimensional conformally flat situation for 
graphene. As said we shall face this problem in detail in the following Section but we anticipate that 
the answer is "yes" . Thus, what is left to do is to identify the metric and define the corresponding 
measuring procedure such that 

p{E, x) = -ImTi (^S'{E, X, x)7°) = -ImTr (^S{E, x, f )7°) = p{E, x) , (3.20) 

and the LDOS for the planar sheet is the same as that for a curved sheet when the two-dimensional 
curvature is such that the 2+1 dimensional metric is conformally flat. 

4 The "nearly two-dimensional" Ansatz for the metric 

The two-dimensional sheet in immersed in a three-dimensional space {x,y,z). One can easily set 
the frame so that the profile of the sheet is z{x, y). Including time then we have the extrinsic frame 
/?" = {t,x,y,z{x,y)) and the intrinsic frame a^ = {t,x,y), where flat (a = 0,1,2,3) and curved 
(^ = 0, 1, 2) indices need to be noticed. The (standard) induced metric procedure gives for the 2 + 
1 dimensional system the metric 

(3)_ d^d^ 



where the choice of four dimensional flat metric, r/afc, is dictated by the special relativistic-like form 
of the action, we just need to remember that the limiting speed for the electronic system is vp and 
not c, all the rest goes through. With this choice for (3°" and a^ the electrons on graphene see a 
spacetime metric of the form 

where a,/3 = 1,2 are the spatial two-dimensional indices and g\i = —1 — (dz/dx)"^, g^-l = —1 — 
{dz/dy)"^, §12 = ^21 ~ —{dz/dx){dz/dy). This we call the "nearly two-dimensional" Ansatz 

(3) 

because, although the metric gliJ is indeed three-dimensional, the curvature content appears to be 

(3) 

all in the two-dimensional spatial part. To see it one first notices that the scalar curvatures for gliJ 
and g^l are equal 

i?(3) = Rm ^ (4.3) 

and then that the Ricci tensors for the two metrics are such that 



^^^^ = J) ^ ■ (4-4) 




These results are due to the Ansatz (j4.2p . On the other hand in any three-dimensional space the 
Riemann tensor (the tensor that bears the full information on curvature) W^' ^aiiv is proportional 
to i?(3) and i2(3)^^ as follows (see, e.g., [20], [21], [22]) 

^(3) f"^^^ ^ eP^\,,^{R^% - \51R^'^) , (4.5) 
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Figure 2: An intrinsically flat bending corresponding to the profile z = cos(3a;). This profile gives 
raise to a flat 2+1 dimensional metric in the appropriate coordinates. 



while in any two-dimensional space (see the Appendix) 



^(2) _ 1 (2)^(2) 



(4.6) 



Thus, from the results (|4.3p - (j4.6p we see that the whole information on the curvature of the three- 
dimensional spacetime with metric ()4.2p is actually in the two-dimensional scalar curvature i?^^', i.e. 
the Riemann tensor only has one independent component, a fact that characterizes two-dimensions, 
for us the {x,y) space of the sheet. This, apparently, means that everything related to curvature is 
happening in the two-dimensional spatial part and we have built on purpose such a situation to be 
in the most conservative position: (i) no curvature on the time part, (ii) the metrics only depend 
on X and y. 

(3) 

Our goal is to see whether a conformally flat g]iv makes sense. Namely, we want to see whether by 
imposing the conditions for three-dimensional conformal flatness we obtain a reasonable condition 
for the shape of the two-dimensional graphene sheet. If that happens then we could reverse the 
argument and say that for such a shape the three-dimensional metric is indeed conformally flat, 
including the time part of it, hence we would be in a truly general relativistic-like system, although 
conformally trivial. What could happen (and we were expecting to happen!) is that only flat sheets 
would be allowed. This case is also of interest, as we shall show, because we are focusing on the 
intrinsic geometric properties, hence, e.g., a shape like that in Fig. [2] is flat. On the other hand if 

(3) 

flat two-dimensional sheets were the only solution then the only conformally flat g\iJ would be the 
flat metric ri\j,y . 

Our problem becomes that of solving the third order partial differential equations 



C, 



^u 



0, 



where 



C^. = eM'=V^i?(^): + e.,.V^i?(3)- 



(4.7) 



(4. 



is (proportional to) the so-called Cotton tensor [2T] [23] and its vanishing is necessary and sufficient 
condition for a three-dimensional metric to be conformally flat. In four dimensions this would not 
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Figure 3: An "egg-tray" shaped (intrinsically curved) sheet corresponding to the profile (I4.9p with 
C = 1 and n = m = 2. We cannot tell whether there is a choice of C,n,m such that this kind of 
profiles give raise to a conformally flat 2+1 dimensional metric. 



be true as it is the vanishing of the Weyl tensor that gives the necessary and sufficient condition 
for the conformal flatness. The latter tensor is identically zero in three dimensions, and this is why 
the Riemann tensor takes the form (|4.5p . 

The task to solve the partial differential equations (j4.7p may become very difficult, hence we 
turned to the help of the package for Mathematica in [21] . With that we have tried with the general 
profile 

z{x,y) = C cos{nx) cos{my) , (4.9) 

with n, ra integers (see Fig. [3]) that could well simulate the ripples found on suspended graphene |14j . 
in the hope to find some values of the amplitudes C and/or of the frequencies n, m for which (14. 7p 
has solutions, besides the fiat solutions n = 0,m = any and n = any, m = 0. The task revealed 
itself to be beyond the capabilities of the program and/or of the machine. Thus we resorted to 
approach the problem using the general argument that, for any given profile z{x,y), there exists a 
coordinate system on the sheet, say it {x, y), such that g^I = —e~'^'^''^'^'6ai3- To find such coordinate 



transformation is surely possible in principle but could be tricky in practice (see, e.g., p5l for some 
examples in a 2+1 dimensional case), but the differential equations (14. 7p surely become much more 
tractable. First, in the coordinates (x, y) the metric (14. 2p becomes 



9l^Jix,y) 



1 











_g-2(T(i,y) 











_^-2a{x,y) 



(4.10) 



With this equations ()4.7p become 

diAa + 2{dia)Aa = (4.11) 

dyAa + 2{dya)Aa = 0, (4.12) 

where A = 9| + 9? is the fiat space two-dimensional Laplace operator. Now we compute R^'^' in 
these coordinates and find (see the Appendix) 

i?(2) = 2e^''Aa . (4.13) 
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By using (14.130 in ()4.1ip and (I4.12p we obtain that the request for gliu in (I4.10p to be conformally 
flat amounts to have a two-dimensional sheet with constant intrinsic curvature R^'^K This, in terms 
of the two-dimensional conformal factor a, amounts to have 



Act = Ke'^" , (4.14) 

with constant Gaussian curvature K = w"^' /2. We wrote the condition on the a in this fashion 
because this way it is easy to recognize it as a classic equation of the geometry of two-dimensional 
surfaces [26j, [27,], namely Liouville's equation [28]. Among the advantages of having such a classic 
equation as our constraint is that all solutions are known [28j . see also [29]-|31). 

Thus we have found that for any graphene sheet that is truly curved, but with constant K, the 
electrons see a truly curved, hut conformally flat, three-dimensional spacetime. To explicitly see it 
one needs to find a third coordinate frame, now involving time (as in principle we cannot and should 
not exclude it) such that 

/ 1 \ 

(4.15) 







' 1 








,lfJ{T,X,Y)-- 


^ g-2S(T,X,y) 





-1 









^ 





-1 



with T, = a only in the trivial case of constant a (and T = e~"t, X = x,Y = y). 

When the two-dimensional spatial metric is such that the Liouville equation ()4.14p is satisfied, we 
have the "most gentle" general relativistic-like situation we could think of, namely the conformally 
trivial one. So, while on the one hand, this means that many characteristics of the flat case 
are untouched in the conformally trivial case (as we explained before), on the other hand some 
features will be general relativistic-like and this, in turn, opens the doors to speculate whether more 
general settings of this kind are meaningful for the physics of electrons on graphene. On the latter 
("exotic") scenario we will briefly comment in the last Section. Let us now look more closely at the 
rich landscape that the constraint (j4.14p opens-up. 

When i^ = we are in the particular case of flat sheets, that clearly gives flat three-dimensional 
metrics. In this case the as permitted are all the infinite solutions of the Laplace equation 

Ad = , (4.16) 

that, as well known, are the harmonic functions in two dimensions. It is also well known that 
in the complex domain, w = x + iy and a{w) = aji{x,y) + ia[{x,y), when an and aj obey the 
Cauchy-Riemann conditions for analyticity, OxCr = dyaj and dxCj = —dyapt, then they are also 
harmonic functions. This means that that are infinite ways we can bend the graphene sheet by 
keeping it intrinsically flat, the representative of one possible family is given in Fig. [2l and that 
this infinite family has tight links with the infinite dimensional Virasoro algebra, the algebra of the 
conformal group in two dimensions, with generators L„ = —w^~^^dw, even if we are dealing with 
a three-dimensional spacetime. On these matters see, e.g., [32]. Of course, to bend the graphene 
sheet it costs elastic energy even though the shape reached is intrinsically flat. Thus, from the 
point of view of the elastic energy, there is no such a thing as a family of invariant shapes. What is 
invariant instead is the way the electrons behave within the intrinsically flat sheets: for them there 
is no distinction among members of the class. 

If we are able to find the coordinate transformation (t, x, y) — )• (T, X, Y) that gives the relation 
between the a of the two-dimensional metric (j4.10p and the T, of the three-dimensional metric 
(I4.15P we would explicitly see how the Virasoro algebra, that surely is there for a, enters the three- 
dimensions via S. This is an interesting mathematical problem on its own and deserves more studies 
[33]. 
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Another interesting point is relative to true curvature for a honeycomb lattice structure as the 
one we have in mind here. As well known, see, e.g., [5], [6], [33] and also [35], for such a lattice 
the curvature can only enter via the appearance of a defect, a fivefold one to have a positive K, a 
sevenfold one to have a negative K. We have not taken this into account in this paper (but shall 
face this in follow-up work [36], [37]) hence, in this respect, we are considering the graphene sheet 
as a continuum. Nonetheless, we want to notice here that if a "no-defect condition" is imposed on 
true graphene then there is no way out to have flat sheets and this, in turn, is not as a limiting 
situation as one might think because the number of allowed shapes is still infinite (and related to 
the Virasoro algebra). 

Thus, within our approximation of curvature not due to defects, the next case to consider is the 
constant K. If we once more move to the complex domain equation (I4.14p becomes 

K 

5^(9iDln(/? = -— 9? , (4.17) 

with 4:dyjdw = d1 + d"^ and ip = e^^'^ . The general solutions are [28] 






^H = 17TrhT77^ ' (4-18) 



with f{w) any analytic function, hence, again, an infinite family of solutions like for the harmonic 
functions. In this case as well this typically two-dimensional rich structure is enjoyed by a three- 
dimensional system but to clearly see how this is reflected on the three-dimensional S one needs to 
find the coordinate transformations [33j {t,x,y) — ?■ {T,X,Y). Nonetheless, this intriguing structure 
underlays our problem and should be explored carefully in relation to possible applications to 
graphene [38]. For instance, vortices have been studied in the context of graphene from various 
perspectives, see, e.g., [39], [40], and we can add to those a new one by considering that among 
the solutions ()4.18p there are (non-topological) vortex solutions [31] given by f{w) = w~'^, with n 
positive integer to which corresponds 



4^2 ^-2(n+l) 

IT (1 + r- 



V'H = — n . .-2n,2 ' (4-19) 



where r = \w\. Current work is devoted to find the actual profile z{x,y) corresponding to those 
solutions and the impact on the density of states |38j. 

5 Conclusions and Discussion 

We can then conclude that, within the approximations explained throughout the papeio, when the 
graphene sheet has a constant Gaussian curvature the three-dimensional spacetime the electrons 
move in is conformally flat and truly curved in all directions, including time, i.e. they experience 
a truly general relativistic environment. Nonetheless, due to the local Weyl symmetry enjoyed by 
the massless Dirac theory, these effects cannot be seen classically nor quantum-mechanically, unless 
"non-invariant measurements" are performed. Geometrically, by "invariant measurements" we mean 
processes of measurement for which times and lengths are given by ds'^ = e~ {dt'^ — dx'^), hence the 
significant vacuum expectations are of the kind (0|O|0) and, due to the invariance, they are equal 
to the fiat vacuum expectations '(0|O'|0)', where O is a generic operator and we use the notation 



^^Besides the low energy limit and the HF-like special relativistic behavior, one crucial approximation here is that 
curvature is not introduced via disclination defects, hence, in this respect we are treating the lattice as a continuum. 
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of the last Section as for conformal factors. Physically, the processes of invariant measurements are 
related to the kind of spacetime S entails. One consequence of the Weyl symmetry leading to this 
general relativistic behaviors is that the electronic LDOS is invariant under the process of curving 
the graphene sheet, provided, as said, the measurements are of the invariant kind and provided 
the two-dimensional curvature is compatible with three-dimensional conformal flatness. For these 
reasons our conclusions in this respect are not necessarily contradicting the results of [3] (see also 
the more recent [41j). 

Of course, this paper is only a preliminary step towards a match with experiments. First, while 
the flat spacetime relativistic-like description of ()3.17p is fully justified from the experimental point 
of view, the same cannot be said for the curved spacetime description of (I3.18P and, in general, this 
is still an open problem. Nevertheless, we have proved here that such curved spacetime description 
indeed makes sense at least for the (infinite) classes of metrics considered here (which physically 
correspond to intrinsically flat and constantly intrinsically curved graphene sheets). For those 
metrics, if one accepts the flat spacetime description then one must accept also a curved spacetime 
description, simply because, through a Weyl redeflnition of the fields, the action is the same. This 
was the main motivation behind our trial metric (14. 2p in the first place: which metric naturally 
gives a curved spacetime graphene action that reduces to the experimentally sound fiat spacetime 
graphene action? 

Second, we need to find the coordinate transformations that give S in terms of the two dimen- 
sional a. Current work is dedicated to this task [38]. There we shall, on the one hand, focus on real 
graphene configurations that can give raise to a three dimensional conformally fiat spacetime, on 
the other hand, we shall be able to appreciate what kind of spacetime is encoded by such bendings 
(besides appreciating the lifting to three dimensions of the Virasoro and Liouville structures [33]). 
In [38] we shall also propose dedicated experiments (or at least the use of existing data) to prove 
(or disprove) at once that the curved space-time action ()3.18p indeed describes electronic transport 
on curved graphene sheet and that Weyl-symmetry is truly in place. It is worth noticing here that, 
although it will surely be very interesting to see how time is affected in the inner reference frame 
through Qqq = exp(— 2S), we have already proved here that it is only a matter of a coordinate change 
(not a change of reference frame) to rewrite the metric where ^qq = 1, i.e. (14. 2p . as exp(— 2S)7/^;y. 
The latter coordinate change always exists if the graphene sheet has constant curvature, and indeed 
it seems a non- impossible task for the experimentalists to realize such sheets (see, e.g., [42] for an 
old overview on generic nanomaterials, and [43j for a more recent work on graphene). 

Once the relevance of conformally fiat spacetimes is established, a natural question is whether it 
makes sense to reverse the logic and consider more general cases of three-dimensional conformally fiat 
spacetimes than those considered here to see which sort of behavior for two-dimensional graphene 
they correspond to ^6J. Intriguing examples are the celebrated conformally fiat black hole in 2-1-1 
dimensions found by Banados, Teitelboim and Zanelli [H] (BTZ) and the kink spacetime of [23] . 
The first case would represent a black hole in a lab [36j, the second case is worth studying on 
its own right but also to explore whether fractionalization of fermion number could occur in such 
topologically non-trivial background ^45j, ^46j, [47j, [36j . 

The latter examples of three-dimensional conformally fiat spacetimes are just two among many 
solutions of the full theory of gravity in three dimensions that was proposed in the early eighties 
in [21] and is still under intense investigation (see [38] for a recent paper). For that theory the 
Euler-Lagrange equations are 

R^, - 2Ag^, + -C^, = , (5.1) 

where all tensors refer to three dimensions and have been defined in the paper (note that this C^i, 
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differs from the one we used in (j4.8p only by a factor proportional to 1/-^, |21j . |23j). A is the 
cosmological constant, and /i is a parameter with the dimensions of a mass. While the first two 
terms on the left side of (15.ip descend from a standard Einstein-Hilbert action (hence are the same 
in any dimension) , the third term descends from an action that is nonzero only in three dimensiona^i 
and it is called the Chern-Simons gravitational term |23j 

Acsig,.) = ^J d'xe>^^' (^r^.a.n, + ^r^^r^^n,) , (5.2) 

where the Christoffel connection can be expressed in terms of the metric, then variations of ^cs 
with respect to g^j^ give the Cotton tensor. Hence, conformally flat spacetimes (C^;/ = 0) are also 
solutions of the full theory {Rfj^u = 2Kg^u). For instance, the kink spacetime mentioned above is a 
solution of Cfj,u = 0, hence of the full theory, while the BTZ black hole is a solution of R^y = 2Ag^y, 
hence, again, of the full theory and conformally flat. Thus, if we intend to hunt for the Weyl invariant 
configurations of graphene depicted here, we should restrict our interest to this class of solutions 
of (j5.ip and answer, case by case, the questions about shapes, invariant measuring procedures, 
invariance of the density of states, etc. Notice that, in general, C^y = does not give the constraint 
of Eqs. ()4.1ip and ()4.12p . hence also nonconstant Gaussian curvature of the two-dimensional sheet 
could give raise to a 2+1 dimensional conformally flat spacetime, hence to a Weyl invariant density 
of states. 

One crucial issue not considered here is about the constraints on the shapes set by the elastic 
properties of the graphene sheet. For instance, in our model, there is no limitation to the amplitudes 
of the intrinsically flat z{x,y) = Ccos(nx) (see Fig. 2) besides the obvious limitations arising 
from the finite size of the sheet (for a square sheet of side L, Cmax = L/2 with n = 1). These 
inputs can only come from the elastic (and dominant) part of the total energy. One direction to 
investigate, then, is to construct a gravity/elastic theory of the graphene sheet along the lines of 
what has been done in [5] and [6], but introducing a further element in that construction that is 
the elastic counterpart of the Chern-Simons gravitational term ()5.2p [37] . Of course to this end we 
need to revisit our analysis because disclination defects are necessary to have a nonzero curvature. 
Nonetheless, if the defects are few and localized (for instance only 12 five-folded defects are necessary 
to have a sphere) and the graphene sample is large enough, the analysis presented here should still 
apply, although less rigorously. 

Let us close this paper with a speculation on the nature of ripples [Il]-[l6]. Since it costs 
extra energy to create a defect it seems reasonable to ask whether ripples are intrinsically flat. We 
could exactly answer this question by using the techniques and invariance (classical and quantum) 
presented here with no approximations because the spacetime is flat: if ripples are intrinsically flat 
the density of states would be fully invariant, while its non-invariant behavior would be simply due 
to a coordinate choice effect (not to an observer effect) and, furthermore, the profile z{x, y) would 
correspond to a two-dimensional conformal factor a{x,y) that is an harmonic function [38j. 
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^It is actually nonzero when n = 4fe + 3, for A: = 0, 1, 2, . 
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Appendix: Geometric conventions and other formulae 

As usual, in n dimension, 



[V^, Vu]-Vx = R^x^^uVp - rVVp^A , (A.l) 

^Vi, = d^Vu — r^^i^Vx and the Riemann and torsion tensors are 

R'^Xij.u = d[LX''n]\ + ^'' [vcj^" ^l]\ : T''^,y = r''[^;y] , (A. 2) 



respectively. By using T^^ = E^{5°'i,dp^ +u}f/^^)e\, which descends from V^e^ = 0, as explained in 
the paper, one has 

{R^,uTb = [di^ ^v\ + ^{p, ^v\) ^ , (A.3) 

and 

^V = E'a{R,^)\e\ , Rx. = E^,{R,,rA , R = E^,{R,,r'E\ , (A.4) 

and 

T''^u = {\i<iu\+^[^.e,X . (A.5) 



where T^ = E^ T"" . For this paper we take T" = (although we keep track of when we require 
that to hold) which gives 

ujabc = \ {E^E:^{d,ebp - d^Cbu) - E^E^{d,eap - d^Cau) + E^E'^id.e,^ - d^e,,)) , (A.6) 

where Uabc = ^a^t^bc- With this (or by using l\2A\\ in the paper and r'*'^,^ in terms of gi^t,y) we have 
the scaling properties of the spin connection 

A, , ab , , ab/^a^a \ , , ab/^a \ J\a Ji] „v /' A '7^ 

^fj. =^fi (ee^j-o;^ (e ^j = e^ ^e V^^ • (A.Yj 

For n = 3, ujfj,ab = ^abc w^^ and uj^ = e^^uj^ and 



w„ 



d 1 bed 



2 



e'"^'* {e^adbEl^ + e^^daE^ + e^^dbE^) , . (A.8) 



For n = 2, uJaab = (-ab^a (or LOa = -^e'^^Wq, ab) and 

e'^^'Ea'^Eb^ = e^P det E = ^e"^ , eabe\e\ = e^^ det e = V^e^/? , (A.9) 

■\/9 

and d^^ijJa\ = ^i3a^^^dsu}-y which give 

1 2 

-R«/37<5 = Vgeapes^e°'''^'di^'UJa' , RpS = —^gpSf^'^d^UJa , R = — e^^S/jWa , (A.IO) 

y9 yj9 



from which it follows 



i^-/3^, = le"/3/7^ , Rps = \gp5R. (A.U) 



2 ' ''" 2 

From these the scaling properties of Wq, , R and Raj3 are 

Ja = oOc,-^e^paf^ , e-''R\p-Rap = gapnga , e" R' - R = 2Uga , (A.12) 
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where primed objects refer to the metric g'^f^ = e^'^gafs and Ogcr = -y=dci{y/gd°^a). Notice that the 

scahng properties of the Ricci and scalar curvature are consistent with R^p = \gapR (multiply both 
sides of the equation for the Ricci by e~'^g°'"). These scaling properties differ slightly from those 
reported in p] (see, e.g., Eq. (28) there). While for the applications of those results to the case in 
point this discrepancy is relevant, hence we took it into account in our calculations here, it is not 
so for the conclusions of [2] because there the important considerations refer to the flat space limit 
for which, e.g., the last relation in (jA.12p becomes Oa = 0, which is the same constraint obtained 
in [2j and the correct condition for conformal symmetry in n = 2. 

Notice that for ^r^^ = — ^q./sc"^'^ it is iOa = e~^'^eQ,^cr^, hence from the last relation in (jA.lOp 

R = 26^" /\a . (A.13) 

Explicit form of the action in the conformally flat case: 



when 

then 

hence from ()A.6P 

From these and from 



g^u = e-2-(-)r?^, , (A.15) 

e;i = e-'^S';, , E^^ = e'^Sa , ^ff = e""'^ , (A.16) 

co^,^ = 5';,{r^abS'',-7]ac6'^W. (A.17) 



71 — 1 

l^'Jab = -^^Ib , (A.18) 



it is easy to see that the action is 



i I (Tx e-("-i)'^ ^I'^ida - - — ^cJa)^ . (A.19) 



When we use ^ = e^a <^(^)\i/' the action reduces to the flat one i f cPx"^' ^^'. 
Let us prove here that S'^{xi,X2) = e^'^"'^^'^''^'^'^^^^'^ S' {xi, X2) satisfies 

i Wx,S^xi,X2) = -^5''{xi-X2) . (A.20) 

First apply the flat Dirac operator to S"^ 

i ^,,5'^(xi,X2) =i^7Va(xi)5'^(xi,X2) + e('^-i)'^("i)<5'^(xi -X2) , (A.21) 

or 

irida - '^aa){xi)S^{xi,X2) = e("-i)'^(^i)5"(xi - X2) , (A.22) 

where we used i ^xiS'{xi,X2) = 6"{xi — X2). Then use the metric in (JA.ISP and ()A.16P - ()A.18P to 
write (1X22]) as 

i7'^e'^(^'i)(a,-^a,)(xi)5'^(xi,X2) = irE!i{df,-^^a^)ixi)S''{xi,X2) (A.23) 



gna 



("^)<5'^(xi-X2) = ^5"(xi-X2), (A.24) 

V9 



i.e. equation ()A.20p . 
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